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Abstract 

In this correspondence, we propose a tight lower bound to the outage probability of discrete-input 
Nakagami-TO block-fading channels. The approach permits an efficient method for numerical evaluation 
q | of the bound, providing an additional tool for system design. The optimal rate-diversity trade-off for the 

Nakagami-rn block-fading channel is also derived and a tight upper bound is obtained for the optimal 

j> . coding gain constant. 
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I. Introduction 

The block-fading channel [1], [2] is a useful channel model for a class of slowly-varying wire- 
less communication channels. The model is particularly relevant for delay-constraint applications 
where channel usage is restricted to only include a finite number of distinct channel blocks, each 
subject to independent flat fading. Frequency-hopping schemes as encountered in the Global 
System for Mobile Communication (GSM) and the Enhanced Data GSM Environment (EDGE), 

K. D. Nguyen and L. K. Rasmussen are with Institute for Telecommunications Research, University 
of South Australia, SPRI Building, Mawson Lakes Blvd., Mawson Lakes SA 5095, Australia, e-mail: 
dangkhoa . nguy en @ postgrads .unisa.edu.au, lars.rasmussen@unisa.edu.au. 

A. Guillen i Fabregas was with the Institute for Telecommunications Research, University of South Australia, Australia. He 
is now with the Department of Engineering, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, UK, e-mail: 
guillen@ieee.org. 

This work has been presented in part at the 2007 IEEE International Symposium on Information Theory, Nice, France, June 
2007. 

This work has been supported by the Australian Research Council under ARC Grants DP0558861 and RN0459498. 



February 5, 2008 



DRAFT 



respectively, as well as orthogonal frequency division multiplexing (OFDM) as encountered in 
more recently proposed wireless communication systems standards can conveniently be modeled 
as block-fading channels. The simplified model is mathematically tractable, while still capturing 
the essential features of the practical slowly-varying fading channels. 

In a block-fading channel, a codeword spans a finite number of B independent fading blocks. 
As the channel relies on particular realizations of the finite number of independent fading 
coefficients, the channel is non-ergodic and therefore not information stable [3], [4]. It follows 
that the Shannon capacity of this channel is zero since there is an irreducible probability that 
a given transmission rate R is not supported by a particular channel realization [1], [2]. This 
probability is named the information outage probability. For sufficiently large codes, the outage 
probability is the lower bound to the word error rate for any coding schemes. 

Considerable efforts have been dedicated to describing the behavior of the word error proba- 
bility and the outage probability for Rayleigh block-fading channels in the high signal-to-noise 
ratio (SNR) regime. In particular, analysis based on worst-case pairwise error probabilities shows 
that at high SNR the achievable word error probability of codes C of rate R (in bits per channel 
use) constructed over a signal constellation X of size \X\ = 2 M behaves as 

SNR-^oc log SNR v ' 



where 



d B (R) = 1 + 



B ( 1 - - 

M 



(2) 



is the Singleton bound [5], [6], [7]. More recently, it has been shown [8] that the optimal SNR 
exponent 

a (R) = sup lim — — — (3) 

v > c *SNR^oo log SNR 

is actually given by the Singleton bound ©. This establishes the Singleton bound as the optimal 
rate-diversity trade-off for transmission over the Rayleigh block-fading channel with discrete 
signal constellations. 

While these results provide significant insight into code design, the analysis techniques do not 
provide explicit tools for the evaluation of the outage probability; a task which usually requires 
extensive numerical computations. To this end, an upper bound to the outage probability of 
Rayleigh and Rician block-fading channels is proposed in [9], [10]. In this paper, we propose 
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a tight lower bound to the outage probability which can be efficiently evaluated for the general 
Nakagami-m block-fading channel [11]. We show that numerical evaluation of the proposed 
bound is very efficient, resulting in significantly less complex computation as compared to Monte 
Carlo simulation. We also show that the optimal rate-diversity trade-off for the Nakagami-m 
fading case is given by d*(R) = md B (R) for any m > 0, and we obtain an upper bound to the 
achievable coding gain for any coding scheme. 

The remainder of the correspondence is organized as follows. In Section II, the system model 
is described for the Nakagami-m block-fading channel, while Section III defines the outage 
probability of this channel. In Section IV, we detail the proposed lower bound for the outage 
probability, as well as an efficient method for the evaluation of the bound. The asymptotic 
behavior of the outage probability is investigated in Section V, where the rate-diversity trade-off 
is extended to include the Nakagami-m fading statistics. Finally, conclusions are given in Section 
VI, while proofs are collected in the Appendices. 

The following notation is used in the paper. Sets are denoted by calligraphic fonts with 
the complement denoted by superscript c. The exponential equality g(£) = £ d indicates that 
Hindoo l °i 9 Jp = d. The exponential inequalities <, > are similarly defined. is the indicator 

function for event [£] ( |_£J ) denotes the smallest (largest) integer greater (smaller) than £, 
and A r l = {£ G A n |£ > 0}. 

II. System Model 

Consider transmission of codewords of length BL coded symbols over a block-fading channel 
with B blocks. Each block is an additive white Gaussian noise (AWGN) channel of L channel 
uses affected by the same flat fading coefficient. The complex baseband expression for the 
received signal is 

y 6 = v/SNR h b x 6 + z b , b=l,...,B, (4) 

where y b G C L is the received signal in block b, x b G C L is the portion of the codeword 
assigned to block b, and z b is a noise vector with independent, identically distributed (i.i.d.) 
circularly symmetric Gaussian entries ~ A/c(0, 1). We define h = (hi, . . . , h B ) G C B as the 
vector of fading coefficients. The fading coefficients are assumed i.i.d. from block to block and 
from codeword to codeword, as well as being perfectly known to the receiver. 
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We consider a channel with a discrete input constellation set X C C of cardinality 2 . Without 
loss of generality, we assume that E[|x| 2 ] = 1, where x £ X, and that the fading coefficients 
are normalized such that E[|/if,| 2 ] = 1. It follows that SNR is the average signal-to-noise ratio at 
the receiver end. Define % = |/i 6 | 2 as the fading power gain. Then, the instantaneous received 
signal-to-noise ratio at block b is 7(,SNR. 

We consider the case where the fading coefficients follow the general Nakagami-m distribution 



[11], [12]. The probability density function (pdf) of \h 



tJ3 



where T(a) is the Gamma function T(a) = J °° t a_1 e~ l dt. It follows that the fading power gain 
76 has the following pdf 

f !! ^rj_ e - m€ e > 

/»(0 = < rM (6) 
I 0, otherwise, 

and cumulative distribution function (cdf) 

F 7b (0 = <^ r(m) - (7) 

I 0, otherwise, 

where T(a,^) is the upper incomplete Gamma function r(a, £) = t a_1 e~ l dt. 

The Nakagami-m distribution represents a large class of fading statistics, including Rayleigh 
fading (by setting m — 1). The distribution also approximates Rician fading with parameter K 
(by setting m = (K + l) 2 /(2i^ + 1)) [12]. Therefore, the proposed analysis for systems with 
Nakagami-m fading is a generalization of previous results in the literature. 

III. Mutual Information and Outage Probability 

The instantaneous input-output mutual information of the block-fading channel with a given 
channel realization h can be expressed as [1] 

B 

B 



1 B 

/(SNR, h) = - J2 

-^AWGN 

(7&SNR), 



6=1 

'Since the complex coefficients ht are perfectly known to the receiver, we can assume phase coherent detection, and thus, 
only the amplitude is affected by the fading statistics. 
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where /awgn(p) is the input-output mutual information of an AWGN channel with SNR p. 
/(SNR, h) is the input-output mutual information of a set of B non-interfering parallel channels, 
each of which is used only for a fraction ^ of the time. When the input signal set X is discrete, 
the mutual information /awgn(p) is given by 



Wn(p) = M-2- a/ ^E 



iog 2 



e 

x'ex 



(8) 



where the expectation over Z ~ A/c(0, 1) can be efficiently computed using the Gauss-Hermite 
quadrature rules [13]. 

Transmission at rate R over the channel in © is considered to be in outage whenever 

1 3 

^^/AWG N (7feSNR) <R. 

b=l 

The corresponding outage probability is given by 

P out (SNR, R) = Pr (± Wn(7&SNR) < r\ . (9) 

IV. Lower Bound to the Outage Probability 

In general, when the channel has a discrete input constellation, evaluation of the outage proba- 
bility in © is complicated since a closed form expression for /awgn(p) is not known. Typically, 
-Pout (SNR, R) is instead evaluated through Monte Carlo simulationsLi which are computationally 
demanding for high SNR. In this section, we propose a lower bound to the outage probability 
with discrete inputs, which can be efficiently computed for any SNR. 

The maximum input-output mutual information for a channel with input signal constellation 
X of size | X | = 2 M is always upper bounded by M. Furthermore, the input-output mutual 
information of the channel can also be upper bounded by that of the channel with Gaussian 

2 Even if the inputs to the channel are Gaussian, for which Jawgn(7(,SNR) = log 2 (l +7bSNR), a closed form expression 
for the outage probability is not known. 
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input. Therefore, for any realization of 7, Jawgn(7&SNR), b = 1, . . . , B is upper bounded by^ 

^awgn(T6SNR) 4 min{M,log 2 (l + 7 6 SNR)} (10) 

log 2 (l + 7b SNR), 7b <w 
M, otherwise 



log 2 (l + 7b SNR), beS c 
M, beS, 



(11) 



where S = |& G {1, 2, . . . , 5} : 7^ > 2 g^~ 1 1 and 5 C denotes its complement. 

Let |<S| be the cardinality of S. Since 7^, b = 1, . . . , B, are independent random variables 
\S\ is a binomially distributed random variable with success rate p = Pr ^ > 2 ^.~ 1 ^ . Hence, 

Pr(|5|=t)= {^jp\l-p) B -\ t = l,2,...,B, 



(12) 



where 

'2 A/ - 1 



p = 1-R 



it 



SNR 



r I m, m 



' SNR 



1 . , (13) 

1 {171) 

Using the upper bound of mutual information in (flOl) and ([111) . we lower bound P out (SNR, R) 
as 

Pl(SNR, R) 4 Pr f I J2 /Xwgn(7.SNR) < flj (14) 



= Pr (^/XwGN(7 b SNR) + ^/XwGN(7 fe SNR) < BJ?) (15) 
\bes bes c J 

= Pr j \S\M + log 2 (l + 7 6 SNR) < Br) . (16) 

Since 75, 6 = 1, . . . , S are i.i.d. random variables, ^ bgi5 c log 2 (l + 7&SNR) is the summation of 
|<S C | = B — \S\ i.i.d. random variables. Each random variable inside the summation is given by 
log 2 (l + 7&SNR) conditioned on b G S c , or equivalently on the event 8, where £ is defined as 



e -^ : -» s w^- (17) 



2«_i 



3 Superscripts u and ^ will denote upper and lower bounds respectively. 
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Denote A b as the random variable log 2 (l + 7feSNR) conditioned on S. Then, the distribution of 
Ab is given by the following proposition. 

Proposition 1: Assume 7 b is a random variable whose distribution is given by ©. Denote 
A b as the random variable log 2 (l + 7f,SNR) conditioned on the event £ given in (TT7T ). The 
distribution of A b is then given by 

I Hg^V ^log(2) p M 

I 0, otherwise. 
Proof: See Appendix U ■ 

Therefore, denoting Ak, k — 1, . . . , |«S C |, as the P — |<S| independent random variables that follow 

the distribution given in (fT8l) . we can write (fl6l) as 



fl-|5| 

Pit (SNR, P) = Pr I \S\M + A k < B R \ ■ (19) 



fe=i 

By conditioning on |<S|, we can express Pg Ut (SNR, R) as 

'b-\s\ 



\S\ = t \ Pt(\S\ = t) (20) 



Pi t (SNR,P) = ]TPr ^ A fc < PP- \S\M 

t=o y fc=i 

B /B-t \ 

= J^Pr I J^A fc < BR-tM J Pr(|5| = t). (21) 

t=0 \fc=l / 

From the distribution in (fT8b . note that Pr(A fc < 0) = 0. Therefore, for any t such that BR — 
tM < 0, or equivalently for all t > \jrf] , the corresponding probability is zero. Hence, we can 
rewrite (TfTI) as 

R?!- 1 /B-t \ 

Pj ut (SNR,P) = Pr [Y, A * <BR - tM Pr (l 5 l=*)- (22) 



t=0 \fc=l 



If we now define the random variable Y t = J2k=i A k, we can write 

rsfl"i_i 

/P\ 

Pl(SNR,P)= £ F yt (PP-iM)(^J^(l-p) B -', (23) 

where F Yt (£,) is the cdf of Y^. 

Since k = 1, . . . , B — t are independent random variables, the pdf of Y t can be evaluated 
by performing B — t convolutions of fA b (C)- Numerically, this convolution can be efficiently 
computed in the frequency domain using fast Fourier transform (FFT) techniques [14]. With this 
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method, we can efficiently evaluate the cdf of Y t , F Yt (£), and therefore we can also efficiently 
evaluate Po Ut (SNR, R) in (|23l) . The evaluation of (|23l) is significantly faster than evaluating 
-Pout (SNR, R) in © using Monte Carlo simulation techniques. 

Numerical results for Nakagami-m block-fading channels with B = 4, M = 4, m = 0.5 
and m = 2 are given in Figure CD The transmission rates considered are R = 1,2,3 bits per 
channel use, which correspond to Singleton bounds (1b(R) = 4,3,2, respectively. The figure 
shows the simulation and analytical curves of the lower bound to the outage probability of the 
channel based on (fl~4)) and (|23l) . respectively, together with the 16-QAM outage simulation curve 
based on ©. We observe that the analytical curves coincide with the corresponding lower bound 
simulation curves. The analytical curves give a tight lower bound to the 16-QAM outage curve. 
Note that the bound is very tight for the important case of R = 1, which, from the Singleton 
bound expression in (0, is the largest rate that can be achieved with full diversity. Figure |2] 
provides a plot of the outage probability of the same channels as a function of the code rate 
R at SNR = lOdB, illustrating the validity of the bound over a wide range of transmission 
rates. Further simulations show that these observations are valid for a wide range of channel 
parameters. We also observe from Figure Q] that the slope of each curve is mdsiR), representing 
the SNR exponent of the outage probability. In the following section, we rigorously prove that 
the optimal SNR-exponent over the channel is 



In proving this result, we characterize not only the SNR-exponent but also the asymptotic coding 
gain. 

V. Asymptotic Behavior 

Using (|23l) and the analysis techniques from [8], we obtain the following result for the 
asymptotic diversity of Nakagami-m block-fading channels, for all m > 0. 

Proposition 2: Assume transmission over the block-fading channel as defined in (HJ) with input 
signal constellation size 2 M . Assume further that the fading power gain 7& is a random variable 
whose distribution is given by ©. In this case, the lower bound on P out (SNR, R) given in (T23T) 
can asymptotically be expressed as 



d*(R) = md B (R). 



(24) 



P^^SNR,^) = /QSNR- mdfl(i?) , 



(25) 
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where d B (R) is the Singleton bound given in ©. Furthermore, /Q is a constant independent of 
SNR given by 

^-F ¥B dB(R) (BR-(B-d B (R))M)( u * (m 5wl^ (fl ^ (26) 

where 



B-d B (R)J (mT(m)) d B(R) 



Proof: See Appendix HB ■ 

This proposition not only shows that the SNR exponent of the outage probability is upper 
bounded by mdsiR) but also gives the asymptotic constant /Q of P"; ut (SNR, R). This is indeed 
useful for code design since it gives an upper bound for the coding gain achieved by any coding 
scheme. At the same time, together with the expression of Po Ut (SNR, R) given in (|23T) . it gives 
a more specific characterization of the outage probability, indicating the word error probability 
(or SNR) region where asymptotic analysis is valid. 

The lower bound to the outage probability and the asymptotic term given in (|25T) are illustrated 
in Figure [3] The same set of parameters as in Figure Q] has been chosen, namely B = 4, M = 
4, m = 2 and R = 1, 2, 3. 

So far, we have shown that d*(R) < mdsiR)- To prove the optimality of the SNR-exponent 
mdsiR), we need to prove the achievability result given in the next proposition. 

Proposition 3: Assume transmission with random codes of rate R and block length L(SNR) 
satisfying 

L(SNR) 

sn™oo log(SNR) = (28) 
over a block-fading channel as defined in (@]) with input signal constellation size 2 M . Further 
assume that the fading power gain 7 b is a random variable whose distribution is given by ©. 
In this case, the SNR-exponent is lower bounded by 

R 



ASMlog(2)(l-f), A< 



d (r) (R) > { M ' ' ^ Mlog(2) 

m(d B (R)-l)+min{m,XMlog(2)(B(l-§)-d B (R) + l)}, A > 

(29) 

Proof: See Appendix [Till ■ 
The preceding propositions lead to the following theorem. 

Theorem 1: Assume transmission over a block-fading channel as defined in (HJ) with input 
constellation size 2 M . Further assume that the fading power gain 7 b is a random variable whose 
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distribution is given by ©. In this case, the optimal SNR-exponent is given by 



d*(R) = md B (R) 



(30) 



for all R, M where B (l - 
Proof: See Appendix [IV] 



) is not an integer. 



As remarked in Appendix [IV] Theorem [TJ can be proved using the methods proposed in [8]. 
However, with the proof proposed here, Propositions [2] and [3] provide additional information. 
In particular, Proposition |2] provides an upper bound on the coding gain /Q, and Proposition [3] 
provides an extension for the SNR-exponent of random codes with finite block length in [8] to 
a more general fading distribution. 

The diversity of random codes for block-fading channels with B = 4, M = 4 and m = 2 is 
illustrated in Figure @J Random codes with block length satisfying A = M f™( 2 ) anc ^ ^ = im iog(2) 
are considered, where A is defined in (|28l) . We observe that the SNR-exponent is always upper 
bounded by mds^R)- Except for points of discontinuity of (Lb{R), the upper bound can be 
achieved by increasing A since S- r \R) and mdsiR) will coincide over larger ranges of R. 



In this correspondence, we have proposed a tight lower bound to the outage probability of 
discrete-input block-fading channels with Nakagami-m fading statistics. The lower bound can 
be computed efficiently and is therefore useful for system design and analysis. We show that 
the optimal rate-diversity trade-off for Nakagami-m block-fading channels is given by m times 
the Singleton bound. We also obtain an upper bound for the achievable coding gain, which is 
useful for code design. 



VI. Conclusion 
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Appendix I 
Distribution and Properties of A b 

Proposition^ Assume % is a random variable whose distribution is given by ©. Denote 
Ab the random variable log 2 (l + 7bSNR) conditioned on the event £ described in (fT7T) . The 
distribution of A b , is given by 



u 


SNR J 


2« Iog(2) 


F -y b ( 


' 2 M_!N 


SNR 




v SNR j 




o, 







fA b (0 = { M^) ' + (3i) 

otherwise. 
Proof: The cdf of is given by 

F Ab (£) = Pr(log 2 (l + 7 b SNR) < CIO (32) 

= Pr (log 2 (l+ 7& SNR)<£ 

Applying Bayes' rule, we obtain 



2 M — 

7 ' £ W > ' 



Pr (j b < 



2«-l ^ 2 M -1 



, SNR ' 1° — SNR , 



If £ < M then 2* - 1 < 2 M - 1 and therefore, 



2 ? - 1 2 M - 1 \ / 2 ? - _ 



Otherwise, if £ > M, 



2«-l 2 M -1\ / 2 M -1\ 

Pr 7fc<T^,7&<^^H = PM76<^r (37) 



SNR ' - SNR J V SNR 

SNR 



2 M 

^ ( "Tnvrr7~ ) ■ (38) 



By inserting (l36l) and (1381) into d34l) . we finally have that 



P T6 I SNR 



= < M^J' (39) 
1, otherwise. 

Now differentiate F Ab (£) in ([39]) with respect to f , noting that = /a 6 (0 and ^^(0 

/ 7i ,(£), we obtain (fTHT) . 
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Proposition 4: Assume ^ b is a random variable whose distribution is given by ©. Assume 
Ab is a random variable as defined in Proposition Q] Asymptotically, the distribution of Ab is 
independent of SNR and is given by 

{ m(2 s_i r -i 2 ; log(2) t <M 
Ifczr ' ' (40) 
0, otherwise. 
Proof: From 06]) and Taylor series expansion, we have 

/ > \ m— 1 

f = - - P SNR 

J <b V SNR j Tim) 



m 



m^2? | \i — I 



r(m) 



-SNR~ (m_1) . (41) 



Similarly, from ©, we have 



F 



/2 M - A . _ r [ m > m lmr) 
V SNR J ~ 1 f(m) 



rfm) ( m 



1 /_2 M -1 



| \ SNR 

r(m) 

m m (2 M - l) m 



mr(m) 

Inserting (gH) and (02) into (HB, we obtain (|4Q|) . 



SNR~ m . (42) 
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Appendix II 
Proof of Proposition [2] 

Define Ak as a random variable described by the distribution function /^ (£) given in (|40l) . 
Further define -Fy t (£) as the cdf of Y t = J2k=i ^h- According to Proposition HI /a 6 (0 = /A b (0> 
and therefore, F Yt {i) = Fpt(0- m addition, Taylor expansion of (fT3l) gives 



P 



P 



T(m) 
m-(2« - 1)- 



= 1, 



(43) 



(44) 



mY\m) 

Since the asymptotic expressions for p, 1 — p and i*V t (£) are finite and non-zero, the asymptotic 
behavior of Po Ut (SNR, i?) in (1231) is found by replacing Fy t (C) with Fy (£), and replacing p, 1 — p 
with their corresponding asymptotic value in (|43l) and (1441) . It follows that 



i* t (SNR,i2) = ^ F Yt (BR-tM) 



t=o 



5 



m m (2 M 



mT(rn) 



B-t 



SNR -m(B-t)_ 



(45) 



Since /a (£) is independent of SNR, Fy t (£) is also independent of SNR. Therefore, the term 
with minimum m(£> — £) dominates the expression in ( I43T ). The dominating term corresponds to 

(46) 



and thus 



B-t=l+ 



M 



- h 



B ' 1 4 

M 



(47) 



which is precisely the Singleton bound. Therefore, we write the asymptotic behavior for (|23l) as 

J^ ut (SNR, i2) = JC t SNR~ mdB {R) , (48) 

where 

(m(2 M - l)) md s(«) 



is independent of SNR. 



B 



B-d B (R)J {mY{m)) d B(R) 



(49) 
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Appendix III 
Proof of Proposition [3] 

The proof follows the same lines as in [8] with the generalization of Rayleigh fading statistic 
to Nakagami-m fading statistic. 

Defining the normalized fading gains as in [15] 

log lb , <m 
ab = ~log(SNR)' (50) 

we have the following result. 

Proposition 5: Assume 75 is a random variable with distribution in ©. Assume further that 
ab is a random variable as defined in (l50l) . In this case, the joint distribution of a = (ai, . . . , as) 
has the following asymptotic behavior 

SNR -™Ef =1 ° 6 fte i« 
/<*(<*)=< + (51) 

0, otherwise. 



Proo/- From 7ft = SNR~ Q \ Therefore, the pdf of a b is 



/a> 6 ) = / 76 (SNR"^) 



dah 



m m 



r(m ; 

The joint distribution of the vector a is then 



r(m) 

SNR- mat exp(-mSNR- ab ) log(SNR). (52) 



/a(a) = (^») B SNR--^-exp (-mgSNR-j . (53) 



It can easily be seen that 



]im tog(/q(a)) _ -mEf=i«6, aeR? (54) 

sn™oo log(SNR) I 0, otherwise. 



Therefore, fct(ot) follows the asymptotic behavior in (pTT) . ■ 
Consider random codes of rate R and block length L = L(SNR) over a signal set of size 2 M 
such that 

A= llm (55) 
SNR-oo log(SNR) 
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Assume the codewords of the code are given by X(i),i = 0, . . . , 2 BLR — 1. Following the 
analysis in [8], the average pairwise error probability between X(0) and X(l) for a given 
channel realization h is given by 



13 



P(X(0)-X(l)|h)< ]J(3t, (56) 



=i 



where (3 b is the Bhattacharrya coefficient 

Pb = 2" 2M £ £ exp (-^J b \x - xf) . (57) 

The union bound of the word error probability for a given fading coefficient is obtained by 
summing over the pairwise error probability of 2 BLR — 1 codewords X(i), i = 1, . . . , 2 BLR — 1. 
Noting that % = SNR 1 "" 6 , we obtain 



P e (SNR|h) 

< exp f -BLM log(2) 



i-£__L 



f>g 2 f 1 + 2~ A/ ]T e -il*- Tsnr-^ j 

6=1 V x^x' / 



(58) 



exp(-5LMlog(2)G(SNR,a)). (59) 



Using (l59l) and the fact that P e (SNR|h) < 1, the average error probability is given by 

P e (SNR) < I min{l,exp(- J BLMlog(2)G(SNR,a))}/ CK (a)rfa. (60) 
J a 

Now, from Proposition [51 we have 

P e (SNR)< I SNR" mE ^i ai 'mm{l,exp(-PLAflog(2)G(SNR,Q : ))}rfQ : . (61) 
Noting that 

lim log 2 ( 1 + 2- M Y e -k\*-*'\W-<*>\ = J °> < 1 (62) 
snr— >oo ^ ^ ; \ M, a 6 > 1, 

we can replace G(SNR, a) in (|6TT) by 

G e (a) = 1-^-1^ i{a n > 1 - 6} (63) 

6=1 

for any e > 0. Therefore, by defining 

B e = \a : G? e (a) < o) (64) 
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and E\ as the complement of B e , we can write (|6Tl) as 

P e (SNR)< I SNR" mE « Q ^a 

JB e nRf 

+ / exp (-log(SNR) (mj^a b + BXM log(2)G e (a) | | da. (65) 

By applying the Varadhan's lemma, the SNR-exponents of the first and second term in (|65l) 
are given by 



B ~\ ( B 

I m 

cteBcnRZ f— * I oteB9 



inf I my a b \ and inf < m) a b + BXM\og(2)G £ (ct) 
respectively. Therefore, the SNR-exponent of the word error probability is given by 
d {r) (R) > sup min < inf mVaj , inf ImS^ a b + BXM\og{2)G e (cx 



e >o | aeB.nRf | f— ' I aez? e c , ... , , 



For the first infimum in (1661) . it can be shown that 

B 



inf < m a b > = m(l — e) B ( 1 ) 



(66) 



(67) 



6=1 

The infimum is attained when \B (l — entries of a. are 1 — e, and the other entries are 
zero. 

The second infimum in (|66l) can be rewritten as 

5 AM log(2) ( 1 - — ^ + inf { V ma b - AM log(2) 1 {a b > 1 - e} j . (68) 
V M J aeB £ c nRS 1^ j 

We consider two cases. If < AMlog(2) < m, the infimum in (|68l ) is zero and achieved when 

a = 0. Therefore, the second infimum in (|66l ) is given by 

5AMlog(2) f 1- ^ J . (69) 



(70) 



M , 

If AM log(2) > m, the infimum in (|68l) is given by 



(m(l-e)- AMlog(2)) 

The infimum is attained when [5 (l — ||)J entries of ck are 1 — e, and the other entries are 
zero. Hence, the second infimum in (l66l) is given by 

R 

1 (;»/( I - r) - A 1/ \oofr\) I! | | - - 

By collecting the results, and noting that the supremum in (|661) is attained when e j 0, we obtain 
the lower bound for the SNR-exponent as in (l2~9l) . 
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Appendix IV 
Proof of Theorem CD 

Clearly P out (SNR, J R) > P Q £ ut (SNR, R), and therefore, 

d* B {R) < md B (R) (72) 

follows from Proposition [2l In addition, by letting L(SNR) — > oo, it follows from Proposition 
[3] that the SNR-exponent mdB(R) is achievable using random codes for all R, M such that 
B (l — -||) is not an integer. ■ 
The theorem can also be proved using the SNR-normalized fading coefficients a 6 = — ^snr.) 
introduced in [15]. The proof given in [8] for the Rayleigh fading case (m = 1) shows that the 
asymptotic behavior of the joint pdf of these coefficients is fcti *-) = SNR - Q & anc [ thus 

d*(R) < inf |XJ a bj = dB ^ (73) 

and 

d*{R) > inf \ Vaj I = d B {R), (74) 
62 Iti J 

whenever £> (l — j^) is not an integer, for some suitably defined sets Bi, i5 2 (see [8] for details). 
In Proposition [5l it is shown that for Nakagami-m distributions the asymptotic behavior of the 
joint pdf of these coefficients behaves as /a(ct) = SNR - " 1 ^^ 1 ab . In this case, the constant 
m factors out from the infimums in (T73T) and (1741) and automatically leads to the desired result. 
While this proof is shorter, Proposition [3] provides the extension of the finite block length results 
of [8], which illustrates the impact of m in the random SNR-exponent d^(R). 
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E s /N (dB) 



Fig. 1. Outage probability of Nakagami-m block-fading channels with B — 4, M — 4, m — 0.5 and m — 2. The thick solid 
lines correspond to the lower bound l !23t . thin dashed lines with circles denote the simulation of l !14t and thin dashed lines with 
squares denote the simulation of ((9) with 16-QAM modulation. 
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Fig. 2. Outage probability for the of Nakagami-m block-fading channels with B — 4, M — 4, SNR = lOdB, m = 0.5 
and m = 2. The solid lines correspond to the lower bound ( 123b . The dashed lines denote the simulation of ((9) with 16-QAM 
modulation. 
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E s /N (dB) 



Fig. 3. Outage probability of Nakagami-m block-fading channels with B = 4, M = 4, m = 0.5 and m = 2. The solid lines 
correspond to the lower bound d23t and the dashed lines to its asymptotic expression given in J25fc using JCe in d26l l. 
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Fig. 4. Optimal and random coding SNR-exponent for Nakagami-m block-fading channels with m — 2, B — 4, M — 4. The 
solid line corresponds to md B (R), dashed-dotted line and dashed line denote the random coding exponent with AM log(2) = 2m 
and AMlog(2) = ^ respectively. 
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